Core-valence correlations for atoms with open shells 
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We present an efficient method of inclusion of the core-valence correlations into the configuration 
interaction (CI) calculations. These correlations take place in the core area where the potential 
of external electrons is approximately constant. A constant potential does not change the core 
electron wave functions and Green's functions. Therefore, all operators describing interaction of M 
valence electrons and N — M core electrons (the core part of the Hartree-Fock Hamiltonian y^-^-^ ^ 
the correlation potential Ei(r, r',i?) and the screening of interaction between valence electrons by 
the core electrons E2) may be calculated with all M valence electrons removed. This allows one 
to avoid subtraction diagrams which make accurate inclusion of the core-valence correlations for 
M > 2 prohibitively complicated. Then the CI Hamiltonian for M valence electrons is calculated 
using orbitals in complete V'^ potential (the mean field produced by all electrons); Ei -I- E2 are 
added to the CI Hamiltonian to account for the core-valence correlations. We calculate Ei and E2 
using many-body perturbation theory in which dominating classes of diagrams are included in all 
orders. We use neutral Xe I and all positive ions up to Xe VIII as a testing ground. We found 
that the core electron density for all these systems is practically the same. Therefore, we use the 
same Ei and E2 to build the CI Hamiltonian in all these systems (M = 1,2,3,4,5,6,7,8). Good 
agreement with experiment for energy levels and Lande factors is demonstrated for all cases from 
Xe I to Xe VIII. 

PACS numbers: PACS: 31.25.Eb, 31.25.Jf 



I. INTRODUCTION 



Accurate calculations for many-electron atoms play 
an important role in many advanced topics of modern 
physics. This includes parity and time invariance vio- 
lating phenomena in atoms search for manifestation 
of possible variation of fundamental constants in astro- 
physical data 0, or in present-day laboratory experi- 
ments Q , improving accuracy of atomic clocks [1] , study 
of super-heavy elements (see, e.g. [5]), etc. Calculations 
are needed for planning of experiments and interpreta- 
tion of the results. 

Atoms of the most interest for the listed topics are 
usually found in the second part of the periodic table 
where measurements or observations are more likely to 
give useful information due to strong enhancement of the 
effects caused by interplay between relativistic and many- 
body effects. On the other hand, accurate treatment of 
relativistic and many-body effects represent a big chal- 
lenge for atomic calculations. Not surprisingly, the num- 
ber of methods capable of producing reliable and accu- 
rate results is very limited. The most advanced methods 
have been developed for atoms with one external electron 
above closed shells. For example, most accurate calcu- 
lations of the parity non-conservation in cesium where 
carried out with two most advanced methods. One was 
the correlation potential (CP) method ^6^ combined with 
the all-order perturbation theory in screened Coulomb 
interaction Q and the other was the linearized coupled 
cluster approach With these two methods, energy 
levels and transition amplitudes for alkali atoms can be 
calculated to the accuracy of fraction of percent while 
hyperfine structure and parity violating amplitudes are 



calculated to the accuracy of 0.5- 1% 0,i,i. 

For atoms with more than one external electron in open 
shells the accuracy of calculations is significantly lower. 
For example, the best accuracy achieved for PNC in Tl is 
around 3% (3% in Ref. [13] and 2.5% in Ref. [ll|)- Typ- 
ical accuracy for energies is about 1% or worse. The 
main challenge is the need for accurate treatment of 
both core-valence and valence-valence correlations. The 
most commonly used methods can be divided in sev- 
eral main groups: (a) many-body perturbation theory 
(MBPT) (see, e.g. [l3|). (b) coupled cluster approach 
(CC) (see, e.g. [ij]), (c) configuration interaction (CI) 
(see, e.g. and (d) multi-configuration Dirac-Fock 

method (MCDF)(see, e.q. [3]). There are also combina- 
tions of these basic techniques. 

All of these method have their limitations. For exam- 
ple, CI usually treats correlations between valence elec- 
trons very accurately but core-valence correlations are 
either totally neglected or only small fraction of them is 
included. MBPT can include more core-valence corre- 
lations, but its application to the correlations between 
valence electrons is limited by the fact that these corre- 
lations are often too strong to be treated perturbatively. 
The CC approach includes certain types of core-valence 
and valence- valence correlations in all-orders and in prin- 
ciple can be formulated for any number of valence elec- 
trons. However, the equations are complicated and most 
of practical realization of the method deal with only one 
or two electrons (or an electron and a hole). 

Significant progress can be achieved by combin- 
ing different techniques. In 1996 we developed a 
method which combines the MBPT with the CI method 
(CI-I-MBPT) [l6[. Here, the second-order MBPT was 
used to construct the effective Hamiltonian in the va- 
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lence space which includes the core-valence correlations. 
It differs from the standard CI Hamiltonian by an ex- 
tra correlation operator E which accounts for the core- 
valence correlations. Single-electron part of this operator 
is very similar to the correlation potential used for atoms 
with one external electron [6]. It was demonstrated that 
inclusion of the core- valence correlations lead to a signif- 
icant improvement of the accuracy of calculations (see, 
e.g. [il,|l7|,[il,[ii). Savukov and co-workers ^ devel- 
oped a version of the method which uses the hole-particle 
formalism. They applied the technique for a calculation 
of the electron structure of the noble-gas atoms [10, Hlj . 

The CI-I-MBPT method was successfully used for 
a number of atoms with two or three valence elec- 
trons [H, [13, [13 (or an electron and a hole [M [HI)- 
Its extension to atoms with more electrons in open shells 
meets some difficulties. It turns out that convergence 
of the MBPT varies very much from atom to atom and 
strongly depends on an initial approximation. The core- 
valence correlations are often too large if treated in the 
same fashion as in our original works and their inclusion 
does not improve the results. 

It is widely accepted that the Hartree-Fock potential 
is the best choice as a zero approximation for consequent 
use of the MBPT due to great reduction of the terms 
caused by exact cancellation between the potential and 
electron-electron Coulomb terms. However, for atoms 
with open shells HP procedure is not defined unambigu- 
ously. This is especially true when the CI method is to 
be used. Here we have freedom of how many electrons 
are to be included into the initial HP procedure and how 
many electrons are to be treated as valence electrons in 
the CI calculations. It was found in our previous work 
[11] that for a wide range of atoms the best choice is the 
so-called V^~'^ approximation. These are the atoms in 
which valence electrons form a separate shell, defined by 
the same principal quantum number. For example, the 
ground state configuration of xenon is [Pd]5s^5p^. Its 
eight outermost electrons have n = 5 while all other elec- 
trons have n < 5. This means that eight outermost elec- 
trons should be treated as valence electrons and the ini- 
tial HP procedure should not include them. This greatly 
simplifies the MBPT, improves its convergence and al- 
lows one to include higher-order correlations in the same 
way as it was done for atoms with one external electron. 

The aim of the present work is to develop a solid the- 
oretical background for use of the y^-^'-f potential as a 
starting point. In principle, this starting point is equiv- 
alent to any other choice of the initial HF potential. 
Indeed, the actual role of the subtraction diagrams in 
the correlation operators T, is to reduce results obtained 
with any zero approximation to results (see an 

explanation below). However, the technique for ]/^-*^ 
is much simpler (no subtraction diagrams) and allows 
us to sum dominating chains of higher order diagrams 
to all orders (it is practically impossible for other choice 
of zero approximation). This results in a higher accu- 
racy. Another advance of the present work is the use 



of a compact basis for valence states. In our previous 
works m, [H, [H, [11] we used the same basis to cal- 
culate S and to do the CI calculations. The basis was 
formed from the eigenstates of the V^-^ potential. This 
had an advantage of having the same single- and double- 
electron matrix elements for all ions of the same atom. 
Moving from ion to ion was easy, requiring only change 
the number of electrons in the CI calculations. How- 
ever, convergence of the CI calculations rapidly deterio- 
rated with growing number of electrons. When number 
of electrons became as large as eight, saturation of the 
basis was very hard to achieve unless huge computer re- 
sources were used. In present work we demonstrate that 
the basis states for valence electrons don't have to be 
eigenstates of the V^~^'^ potential. Instead, HF states 
calculated in the mean field of all electrons (of a neutral 
atom or corresponding positive ion) can be used after mi- 
nor modifications. In this case we have to recalculate the 
CI basis when we change the number of valence electrons 
M. However, the gain is much larger. Since HF states 
are already good approximations to the wave functions of 
valence electrons we can limit the basis to just few states 
in each partial wave. Therefore, even for eight valence 
electrons the CI matrix is small, its calculation and diag- 
onalization takes little time but the final results are very 
accurate. 

We calculate energy levels and Lande g-factors for neu- 
tral xenon and all its positive ions from X II to Xc VIII 
for illustration on how the technique works. Good agree- 
ment with experiment is demonstrated for all cases while 
very little computer resources are needed on every stage 
of the calculations. 



II. CORE ELECTRON DENSITY AND 
POTENTIAL IN v'^-"' AND 
APPROXIMATIONS 

The effective Hamiltonian of the CI method has the 
form (see, e.g. (llj) 

1=1 i<j ' "* ' 

Summation goes over valence electrons, h\{ri) is the one- 
electron part of the Hamiltonian 

hi = cap + (/3 - l)mc^ + Vcore, (2) 

r 

a and (3 are the Dirac matrices, Ze is the nuclear charge 
and Vcore is the electrostatic potential created by the core 
electrons. Regardless of initial approximation used to 
calculate core and valence states, the valence electrons 
never contribute to Vcore directly. They can only con- 
tribute to Vcore via the self-consistent HP procedure or 
via any other potential used to represent valence elec- 
trons. If the core electrons and valence electrons belong 
to different shells the effect of the valence electrons on 
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n = (r + (S log(r))/h 

FIG. 1: Electron densities (multiplied by r^) of Xe I and 
Xe IX. Atomic core (n=l,2,3,4) of Xe I, solid line; 5s and 5p, 
dashed line; electron density of Xe IX, dots. 



electron states in the core and thus on Vcore can be ex- 
tremely small. Indeed, in this case the overlap between 
density of the valence electrons and density of the core 
electrons is small. Therefore, the exchange interaction 
between the core and valence electrons, which is propor- 
tional to the overlap, is negligible in comparison with en- 
ergy of the core electrons. On the other hand, the direct 
potential created by the valence electrons is practically 
constant inside the core since nearly all charge of the va- 
lence electrons is located outside the core. Constant po- 
tential corresponds to zero electric field and cannot have 
any effect on the wave function of the core electrons. The 
only effect of the constant potential Vq is in energy shift 
6E = Vq. However, it does not change the single-particle 
wave functions and Green's functions of core electrons 
since the wave equation contains the difference E — Vo 
which does not change. We may formulate this conclu- 
sion using the perturbation theory. In first order in Vq 
a core state a in the approximation and a in the 

approximation are related by 



(3) 



If the potential Vq is constant, the matrix element 
(a|Vb|n) = Vb(a|n) = due to the orthogonality con- 
dition. This explains why the changes of the core wave 
functions, density and potential are very small. 

Small overlap between the core and valence states usu- 
ally takes place when these states correspond to different 
atomic shells defined by the principal quantum number 
(see, e.g. [H,[23,|2^). In case of xenon, eight outermost 
electrons have principal quantum number n = 5 (the 



5s^5p^ ground state configuration), while all core elec- 
trons have n < 5. Therefore, if the eight electrons are 
considered to be valence electrons we should expect that 
they have little effect on the core states. Fig. [T] shows 
electron densities of Xe I and Xe IX calculated in the 
HF approximation. For the neutral Xe I electron densi- 
ties of valence and core electrons are shown separately. 
One can see that the overlap between them is indeed very 
small. Therefore, it turns out that when electron density 
of Xe IX is calculated it practically coincides with the 
electron density of the core states of neutral Xe. The 
former is shown by dots of Fig. [T] Resolution of this fig- 
ure doesn't allow us to see any difference between elec- 
tron densities of Xe IX and the core of Xe I. This is in 
spite of huge difference in energies of core states of two 
atoms. One may argue that huge difference in energies 
should lead to a noticeable difference in wave functions, 
at least on large distances. Indeed, a wave function of an 
atomic electron has asymptotic defined by its energy and 
potential 
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(4) 



However, in the area up to the radius of the valence shell 
there is actually no difference in _E — V(r) for the core 
orbitals in Xe IX and Xe I since 8E = l^W). The differ- 
ence in asymptotic behavior appears only near the radius 
of the valence shell where the core electron density is ex- 
tremely small. 

Thus we conclude that the core electron density and 
potential have practically no dependence on the number 
of valence electrons if the valence electrons are in a dif- 
ferent shell. 



III. THE CORE- VALENCE CORRELATION 
CORRECTIONS IN THE AND 1/^ *^ 
APPROXIMATIONS 

The use of the Feynman diagram technique allows us to 
express the core- valence correlation corrections in terms 
of the single-particle wave functions and Green's func- 
tions 7] (see Appendix). Therefore, all the arguments 
presented above are applicable when we consider calcu- 
lation of the correlation operators Si and £2; they may 
be calculated using y^-^ basis for core electrons. 

It may be instructive to clarify this conclusion using 
more popular Schroedinger perturbation theory where 
explicit summation over intermediate states is involved. 
The correlations between the valence and core electrons 
as well as the screening of the interaction between the 
valence electrons happen inside the area occupied by the 
core electrons. Let us enclose the core by a sphere with 
zero boundary condition for the core electrons. This al- 
lows us to reduce the core electron problem to the discrete 
spectrum. Let us now consider the interaction of the core 
electrons with external electrons using the perturbation 
theory. The constant potential Vq of external electrons 
does not change the core electron wave functions. It also 
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does not change the energy differences £"„ — Em, between 
the enclosed "core" states, they are shifted by the same 
energy Vq (note that these enclosed states form complete 
basis set inside the sphere). Therefore, all the terms in 
the perturbation theory for the core-valence interaction 
(beyond the mean field which we take into account in 
the valence orbitals) do not depend on the specta- 
tor valence electrons. This is why we can calculate all 
core-valence correlations using the y^-*^ core orbitals. 
To avoid misunderstanding we should note that we use 
this picture for the explanation only, no special bound- 
ary conditions for core electrons are needed for actual 
calculations (it is obvious if we use the Green's function 
technique; all the integrals over coordinates are domi- 
nated by the core area where the correlations between 
the valence and core electrons actually happen). Note 
that below we do not neglect effects of Vq, we only treat 
them as a perturbation since the non-diagonal matrix el- 
ements (a|Vo|n) are small. 

The effective Hamiltonian of the CI+MBPT method 
has the form similar to ([T]) but with extra terms for sin- 
gle and double electron parts of it. These terms, for 
which we use notation E, describe the core-valence cor- 
relations [iBl- There is a single-electron operator Si 
which is added to the single-electron part hi ([2]) of the 
CI Hamiltonian: 
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hi{r) —f hi+ El. 



(5) 



El describes correlations between a particular valence 
electron and core electrons. It is very similar to the cor- 
relation potential E used for atoms with one external 
electron (see, e.g. Q). 

There is also a two-electron operator E2 which modifies 
Coulomb interaction between valence electrons: 
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FIG. 2: Second-order diagrams for single-electron correlation 
operator 
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Physical interpretation of E2 is the screening of Coulomb 
interaction between valence electrons by core electrons. 

When number of valence electrons is greater than 2 
there is also a three-body operator E3 [l^ and higher- 
order many-body operators S4, E5, etc. . However, they 
are usually very small and we will not consider them in 
the present work. 

The full set of diagrams for Ei and S2 in the sec- 
ond order of MBPT is presented on Figs. [2 [3l [4] and [3 
It contains the so-called subtraction diagrams which are 
proportional to Vcore — V^^ , where Vcore is the potential 
of the core electrons as in the effective CI Hamiltonian 

^ yHF -j-jj^g potential in which states of the core were 
calculated. Note that subtraction diagrams vanish in the 
yN-M approximation: Vcore — V^^ . 

The origin of the subtraction diagrams is clear from 
the definition of the perturbation (residual interaction) 
operator U = Hexact — Hq where H^xact is the exact 
Hamiltonian and Hq is the zero approximation Hamilto- 
nian. If the field of external electrons is included in Hq it 
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FIG. 4: Second-order diagrams for double-electron correlation 
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FIG. 5: Subtraction diagrams for E, 
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produces additional contributions which we call the sub- 
traction diagrams. Thus, the potential Vq appears with 
positive sign in the mean field (and core wave functions), 
and with negative sign in the residual interaction (and 
subtraction diagrams). If we calculate all correlations 
exactly, to all orders, Vq must disappear in the final re- 
sult. In any finite order of the many-body perturbation 
theory there are only partial cancellations; lower orders 
of expansion in Vq are canceled out. Thus, the role of the 
subtraction diagrams is to cancel the potential of specta- 
tor valence electrons acting on the core electrons (effect 
of valence electrons on the core lines of the diagrams). 
In other words, the subtraction diagrams guarantee that 
in anj given order of expansion in Vq the operators Si 
and E2 are reduced to the results of y^^*^ approxima- 
tion. Therefore, if the non-diagonal matrix elements of 
Vq are small the v'^"^^ approximation is the best zero 
approximation since the calculations arc much simpler 
(no subtraction diagrams). 

It is easy to see all these cancellations of Vq explicitly, 
order by order in Vq. Here one should remember that 
change of the valence electron energies due to change of 
the core Hartree-Fock potential (which formally has the 
first order in the Coulomb interaction) is actually can- 
celed by the second order subtraction diagrams; contri- 
bution of Vq into the core wave functions in the second 
order diagrams is canceled by the third order subtraction 
diagrams, etc. 

IV. HIGHER ORDER TERMS IN E 

We have seen above that if the electrostatic potential 
Vq created by valence electrons is nearly constant inside 
the core then the v^~'^^ approximation is equivalent to 
the v^~^ + Vq approximation, where Vq can have con- 
tributions from all M valence electrons, or any fraction 
on them or it can be just a model potential. The only 
condition is that Vq is nearly constant inside the core. 
This means that without any compromise on accuracy 
we can do the calculation in the approximation 
which is technically more simple. Another advantage of 
using v^~^^ is that we have the same effective Hamil- 
tonian for any number of valence electrons from 1 to M. 
Therefore we can do the calculations in a very similar 
way for all corresponding ions as well as for a neutral 
atom. 

Eliminating subtracting diagrams in the v^~^^ ap- 
proximation makes Si practically identical to the cor- 
relation potential S used for atoms with one external 
electron. Therefore, we can try to improve the accu- 
racy of calculations by including important higher order 
terms into Si the same way as it was done in a num- 
ber of calculations for alkah atoms (see, e.g. 0, Q)- We 
include two dominating classes of higher order diagrams 
into calculation of Si. One is screening of the Coulomb 
interaction between valence and core electrons by other 
core electrons. Another is interaction between an elec- 



TABLE I: Removal energies of lowest states of Xe VIII (cm ^) 
in different approximations; comparison with experiment 



State 


HF 


f;{2) 


^{00) 


Expt [26] 


5Si/2 


839764 


858722 


854842 


854755 


5Pl/2 


725342 


741343 


738280 


738288 


5p3/2 


707377 


722378 


719550 


719703 


5^3/2 


536494 


546174 


544092 


544867 


5(^5/2 


533632 


543150 


541096 


541939 


4/5/2 


572050 


590725 


587324 


589594 


4/7/2 


571684 


590056 


586717 


589044 



tron excited from the core and a hole in the core created 
by this excitation. Both classes of diagrams are included 
in all orders (see, e.g. 0,0 for details). 

We use notation S^°°-* for the all-order S operator as 
compared to Y,\ for the second-order operator. The ef- 
fect of inclusion of second and higher order correlations 
can be illustrated by calculating of the energy levels of 
Xe VIII. This ion has only one valence electron and cal- 
culations for it can be done the same way as for other 
single- valence electron atoms (see, e.g. 0]). Instead of di- 
agonalizing the CI matrix we solve HF-like equations for 
valence electrons in coordinate space, with Si included 
in it: 

(/ii + Si - e)V = 0. (7) 

Here single-electron Hamiltonian hi is given by ^ while 
2ji can be either sf^ or If no Si is included then 

eq. ([7]) gives HF energies and wave functions. 

The results of calculations are presented in Table U and 

compared with experiment. One can see systematic sig- 

" (2) 

nificant improvement of the results when first Y.\ and 
then sj""-* are included. 

We are now going to use the same S^""-* operator for 
all ions from Xe VII to Xe II and for neutral xenon. For 
all these ions which have more than one valence elec- 
trons the S2 operator should also be included. In the 
yN-M appj-Qximation the S2 term is given by diagrams 
on Fig. m and no subtraction diagrams are needed. To 
include higher-order correlations into S2 we use screen- 
ing factors the same way as we do this for the exchange 
diagrams of Si (Figs. [12 and[ll4) (see, e.g. 0). To ex- 
plain how screening factors are found and used we need 
to go into more details on how the all-order correlation 
operator S^°°'* is calculated. We use Feynman diagram 
technique to calculate direct diagrams (Fig. d 1 and 3). 
It allows us to include an infinite chain of screening di- 
agrams in all orders [Tj]. Application of the Feynman 
diagram technique to exchange diagrams (Fig.[2]2 and 4) 
is much more complicated Q. On the other hand these 
diagrams are usually an order of magnitude smaller than 
direct diagrams. Therefore it makes sense to use an ap- 
proximate method by introducing screening factors. We 
assume that screening of Coulomb interaction between 
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TABLE II; Basis states of valence electrons used in the CI 
calculations, their total number (A^) for each atom or ion and 
HF configurations in which they were calculated 



Atom 


N 


Basis states 


Confif 


durations 


Xe I 


15 


5s,6s,7s,5p,6p,7p,5d,6d,4f 


5s"5p^ 


5s^5p''n/ 


Xe II 


10 


5s,6s,5p,6p,5d,4f 


5s2 5p^ 




Xe III 


8 


5s,6s,5p,6p,5d 


53^5/, 




Xe IV 


10 


5s,6s,5p,6p,5d,4f 


5s2 5p^ 




Xe V 


8 


5s,6s,5p,6p,5d 


5s25p^ 


5s^5pnl 


Xe VI 


8 


5s,6s,5p,6p,5d 


5s^5p, 




Xe VII 


8 


5s,6s,5p,5d 


5s^ 


5snl 



core and valence electrons depends only on multipolar- 
ity k of Coulomb interaction. Screening factors fk are 
calculated as ratios of partial contributions to Ei: 

A- = (i]i°°V(sf>, (8) 

where only direct diagrams are included in S^'^'' and 
and only screening of Coulomb interaction but no hole 
particle interaction is included in . The values of fk 
found from calculations for alkali atoms are 

fo = 0.72, /i = 0.62, h = 0.83, 

/3 = 0.89, /4 = 0.94, /s = 1.00, etc. 

The values of fk change very little from atom to atom 
and the values presented above can be used for xenon. 
This is supported by the results obtained for Xe VIII 
(see Table HI). 

The effect of S2 on atomic energies is much smaller 
than those of Ei. Therefore we can also treat higher- 
order correlations in E2 in an approximate way, via 
screening factors, as we do this for exchange part of Ei. 
We replace every Coulomb integral Qk on all diagrams 
on Fig. m except diagram Fig. U)! by its screened values 
fkQk where screened factors fk are taken as in For 
the diagram Fig. [Hi only one of the Coulomb integrals 
is replaced by its screened value. This is because this di- 
agram can generate only one infinite chain of loops rep- 
resenting screening. Therefore, screening should by in- 
cluded only once. This this very similar to the all-order 
treatment of the direct diagram for Ei (Fig.[2ll and[2l3). 
If this diagram is expressed in terms of screened Coulomb 
interaction, only one of two Coulomb integrals should be 
replaced by a screened one (see [Tll^] for details). 

V. BASIS 

There are two single-electron basis sets in this prob- 
lem. One is used to calculate S and other is used to 
construct many-electron states of valence electrons for 
the CI calculations. 

In principle, it is possible to use the same basis for both 
purposes and we did so in many of our earlier calcula- 
tions [13, mil mill. The most convenient choice for 



the basis is the basis consisting of single-electron states 
calculated in the V^~^^ potential. We use B-spline tech- 
nique to calculate the basis. Lower and upper compo- 
nent of each basis set is expressed as linear combination 
of 40 B-splines in the cavity of radius of 40aB. Expan- 
sion coefficients are found from the condition that the 
basis states are the eigenstates of HF Hamiltonian ([2]) 
with the V'^~^^ potential. The advantages of this ap- 
proach are many: core and valence states are orthogonal 
automatically, the basis is reasonably complete and does 
not depend on number of valence electrons. The latter 
means in particular that if we want to change number 
of valence electrons (e.g. to do calculation for another 
ion) we don't have to recalculate single and double elec- 
tron matrix elements. The shortcoming of this approach 
is rapid increase of the size of the CI matrix with the 
number of valence electrons. Indeed, typical number of 
single-electrons basis states needed to get saturation of 
the basis is around 100. The number of ways, valence 
electrons can be distributed over these 100 states grows 
very fast with the number of valence electrons. For eight 
electrons like for Xe I the matrix reaches unmanageable 
size, even when some configuration selection technique is 
used. 

In present work we use the basis described above only 
for calculation of E. For the CI calculations we use very 
compact basis of HF states of corresponding ion or neu- 
tral atom. For example, we perform HF calculations for 
neutral Xe I in its ground state [Pb]5s^5p^ in the V'^ 
approximation and then use the 5s and bp states as the 
basis states for the CI calculations for Xe I in the V'^~^ 
approximation. Other basis states like 6s, 6p, etc. are 
obtained by removing one 5p electron from the atom and 
calculating these states in the frozen field of remaining 
electrons. The states obtained this way are not orthog- 
onal to the core which corresponds to the y^^^ poten- 
tial. In the relativistic case the basis states for valence 
electrons must also be orthogonal to the negative energy 
states (positron states). Both conditions (orthogonal- 
ity to the core and to the negative energy states) can 
be achieved by projecting a basis state on the B-spline 
states above the core: 

\v)^\v')=Y,\^){^\v). (10) 

i 

Here summation goes over states above the core. Func- 
tions \v') are more suitable for the CI calculations than 
states because they don't have admixture of the core 
and negative energy states. 

If more than one state of particular symmetry is in- 
cluded into the basis (like, e.g. the 6p and 7p states for 
Xe I) they also need to be orthogonalized to each other. 

Full list of valence states for xenon and its ions used in 
the calculations are presented in Table [TTl First column 
shows an atom or ion, second column gives total number 
of valence basis states, then states are listed together 
with the configurations in which they were calculated. 
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Note that every state with I > consists of two function, 
e.g. 6p stands for 6pi/2 and 6^3/2, etc. Note also that the 
number of basis states is always small, much smaller than 
about 100 needed with the B-spline basis. This greatly 
overweights an inconvenience of recalculating the basis 
for every ion or atom. 



VI. CALCULATIONS FOR XENON AND ITS 
IONS 

In this section we present calculations for xenon and 
its ions. The whole calculation scheme consists of the 
following steps (we use Xe I as an example): 

1. HF for Xe IX, potential is obtained. 

2. Calculation of B-spline states in the V'^~^ poten- 
tial. 

3. Calculation of Si. 

4. HF for Xe I, the 5s and 5p basis states are obtained. 

5. Calculation of valence basis states. 

6. Calculation of single and double-electron matrix el- 
ements, including matrix elements of S2. 

7. Calculation and diagonalization of the CI matrix. 

Ar first glance this scheme doesn't look very simple. 
However, none of the steps listed above are very time 
consuming or require large computer power. The most 
time consuming step is calculation of S (Ei in step 3 and 
E2 in step 6). An efficient way of calculating both Si and 
S2 is presented in the appendix. The timescale to obtain 
all results presented in this section while using a PC or 
a laptop is one day. 

Results for neutral xenon are presented in Table IIIII 
while results for six positive ions from Xe VII to Xe II 
are presented in Tables [iVl El IVH IVIIIVIIII and [iXl For 
neutral xenon (Table IIIip we study in detail the role of 
core-valence correlations by including them in different 
approximations. The basis for valence states is kept the 
same in all cases (see previous section for the description 
of the basis) . The approximations are 

1. First, we present the results of the standard CI 
method, with no core-valence correlations (the 
"CI" column of Table IIIip . Accuracy for the en- 
ergies as compared to experimental values are not 
very good. However the difference does not exceed 
10% which is sufficiently good for many applica- 
tions. This is in spite of the fact that calculations 
for neutral xenon were done with atomic core cor- 
responding to highly ionized Xe IX. This is another 
confirmation that change in the core potential Vcore 
from Xe I to Xe IX is very small. 



^ (2) 

2. Second-order Y,\ is added to the effective Hamilto- 
nian (the "sf^" column of Table HIJ). The resuhs 
are significantly closer to the experiment but the 
correction is too large. This is similar to what usu- 
ally takes place with the second-order correlation 
correction for atoms with one external electron. 

3. Second-order t,^^^ is also added (the "E^^feS^^^" 

" (2) 

column). As one can see ^^ts in opposite di- 

^ (2) 

rection to Ei and the results are even closer to 
the experiment. 

4. Higher orders are included in Ei while E2 is not 
included at all (the column). The effect 
of higher orders in Ei is numerically close to the 
effect of E2 as is evident from the comparison with 
previous column. This coincidence is accidental. 

5. Higher orders are included in Ei while E2 is in- 
cluded in second order (the "eJ°°''&E2^''" column). 
The results are improved but for many states the 
correction is too large. 

6. Higher orders are included in both Ei and E2 (the 

column). This is the most complete calcu- 
lation we have in present work. Here we also in- 
cluded calculated values of Lande's (7-factors. The 
g-factors are very useful for identification of the 
states, especially for atoms with dense spectrum 
where calculations do not always reproduce the cor- 
rect order of the levels. 

The last column of Table IIIII presents the difference be- 
tween experimental and calculated energies where calcu- 
lated energies correspond to the most complete calcula- 
tion (S(°°^): A = Eexpt — Ecaic- This difference does 
not exceed 2% and should mostly be atributed to incom- 
pleteness of the basis. Indeed, it is hard to expect that 
the basis consisting of only 15 single-electron states (from 
one to three in each partial wave from I — to I — 3) 
to be complete. Test calculations show that adding more 
states to the basis do have some effect on the energies 
of the states. The effect is larger for higher states. For 
example, it is hard to expect any reasonable accuracy for 
the states of the 5s^5p^6d configuration without having 
the 6d state in the basis. But adding the 6d state to the 
basis also have some effect on the lower 5s^5p^5c? config- 
uration. The detailed study of the ways to saturate the 
basis goes beyond the scope of the present work. 

Tables llVllVl lVll I VIII VlHl and HXl present our results 
for xenon positive ions from Xe VII to Xe II. Only re- 
sults obtained in the "best" approximation (E(°°)) are 
included. Calculations for the ions start from point 4 
in the scheme presented in the beginning of this section. 
This is because first 3 points are exactly the same as for 
neutral xenon. Note that one of the most time consum- 
ing steps, calculation of Si, doesn't need to be repeated. 
Basis states for valence electrons used in the CI calcula- 
tions are described in previous section (see Table |TT| . We 
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TABLE III: Ground state removal energy (a.u.), excitation energies (cm ^) and ,g-factors of lowest states of Xe in different 
approximations 



State 


J 


CI 








S<°°'&E^^^ 


E' 


;oo) 


Expt 










E 


E 


E 


E 


E 


E 


9 




9 


A 


OS op 







-15.21 


-15.76 


-15.69 


-15.53 


-15.48 


-15.49 




-10. Di 






OS op OS 




2 


62710 


70595 


68587 


68289 


66310 


67040 


1.4994 


( UDo 


i.ouuyo 


28 






1 


64013 


71916 


69873 


69594 


67573 


68319 


1.2157 


DoU40 


i.zUOO 


-274 


OS op OS 







71616 


80192 


78425 


77568 


75824 


76480 





( Diy 1 




-283 






1 


72896 


81586 


79764 


78925 


77120 


77799 


1.3160 


I I loO 


1 Q01 


-614 


2r 5^ 

OS op op 




1 


72707 


81332 


79318 


78519 


76533 


77283 


1.8559 


770t?n 


l.OOZ 


-14 









76219 


84202 


82183 


81602 


79607 


80350 


0.0000 


ouiiy 




-231 


OS Op op 


[5/2J 


2 


73810 


82369 


80304 


79587 


77545 


78307 


1.1005 


7Q 1 on 


1 111 r\Q 
i.iilUo 


-187 






3 


74045 


82627 


80568 


79837 


77802 


78564 


1.3333 


7Q /IflQ 


i.ooD 


-161 


OS op op 


2 Fq /ol 
[6/2\ 


1 


74755 


83339 


81261 


80549 


78496 


79260 


1.0232 


7Qnc^f; 

( oyoD 


i.U2o4o 


-304 






2 


74964 


83540 


81466 


80755 


78705 


79468 


1.3913 


7001 


1 QQOa 

i.oooD 


-256 


OS Op oa 







76068 


84821 


82949 


82082 


80235 


80919 





7n771 

r y M 1 




-1148 






1 


76259 


84946 


83067 


82221 


80367 


81054 


1.3786 


( yyf5 ( 


i.oyO 


-1067 


C «2 r „5 r J 

OS op oa 


2 [rj /Ol O 

[7/2J 


4 


76425 


84622 


82869 


81997 


80257 


80900 


1.2500 


OA1 n7 


i.zoUb 


-703 






3 


77283 


85513 


83764 


82877 


81141 


81783 


1.0762 


ouy ( u 


1 nv/io 
i.u / 4y 


-813 


OS Op oc( 


2 Fq /Olo 


2 


76370 


84667 


82973 


82001 


80323 


80947 


1.3775 


CAQOQ 


L.oiOO 


-624 






1 


80595 


89175 


86989 


86445 


84279 


85087 


0.9900 


QQOQA 

oooyu 




-1197 


K ^2 r „5 r J 

OS op oa 


2 fp: /ol o 
[5/2J 


2 


78310 


86627 


84853 


83959 


82199 


82846 


0.9419 


o 1 nof? 
OiyzD 




-920 






3 


78626 


86996 


85209 


84306 


82531 


83184 


1.2179 


oz4oU 




-754 


OS Op iS 


2 Fq /olo 


2 


80504 


89101 


86967 


86367 


84254 


85042 


1.4910 


ooioy 




147 






1 


81064 


89742 


87648 


86951 


84880 


85655 


0.9759 


o044U 




-215 


C ^2 r „5 '-7„ 

OS Op (p 


2Fi /Ol 

[1/2J 


1 


83048 


91935 


89859 


89051 


87001 


87773 


1.7930 


07n07 


1 '70'70 


154 









84221 


92910 


90834 


90110 


88057 


88825 





QQQ/1 




17 


OS op ip 


[5/2J 


2 


83464 


92252 


90160 


89397 


87330 


88104 


1.1107 


CQQ c;o 

ooooz 


i.i2 / D 


248 






3 


83558 


92350 


90262 


89494 


87429 


88204 


1.3333 


oQ/1 art 


i.ooU 


265 


OS op 6p 




1 


83700 


92521 


90429 


89662 


87594 


88369 


1.0216 


Q Q Q7n 


A 700 p: 


10 






2 


84883 


94445 


92365 


91398 


89352 


90122 


1.1497 


OQ1 f?0 

oy iDZ 


1 1 on 


-960 


5s^5p^6d 


2 [1/2]° 





83637 


92398 


90367 


89562 


87557 


88309 





88491 




182 






1 


83728 


92489 


90466 


89654 


87658 


88405 


1.3430 


88550 




145 




2 [3/2] 


2 


83832 


92610 


90519 


89759 


87693 


88466 


1.3843 


88687 


1.3520 


221 






1 


84392 


94128 


92055 


90956 


88909 


89680 


0.6345 


88745 


0.9039 


-935 




2 [3/2]° 


2 


83947 


92702 


90695 


89870 


87892 


88632 


1.3165 


88708 




76 






1 


86666 


95327 


93442 


92542 


90699 


91378 


0.6980 


90032 




-1346 




2 [7/2]° 


4 


84071 


92767 


90743 


89952 


87956 


88701 


1.2500 


88912 




211 






3 


84273 


92937 


90939 


90134 


88167 


88900 


1.0926 


89025 




125 




2 [5/2]° 


2 


84610 


93238 


91272 


90449 


88514 


89232 


0.9548 


89243 




11 






3 


84900 


93476 


91533 


90700 


88790 


89496 


1.2085 


89535 




39 



VII. SOME SPECIAL CASES 

A. Highly excited states 

One of the additional advantages of the use of ba- 
sis for valence states is the possibility to study highly 
excited states with a very short basis. To get to a highly 
excited state with an universal basis like B-splines one 
has to calculate all states of the same parity and total 
momentum J which are below of the state of interest. 
Also, the completeness of the basis deteriorates rapidly 
while going higher in the spectrum. basis is free 

from these problems. To calculate highly excited states 
of a particular configuration it is sufficient to include into 
single-electron basis for valence states only states which 
correspond to this configuration. For example, the states 



use shorter basis for the ions because we calculate only 
lowest states. To go up in the spectrum we would need 
to extend the basis similar to what is done for Xe I. The 
analysis of the data in TablesHVl FVl FVll IVIIIVHH andHXl 
show that the accuracy is generally very good in spite of 
very short basis. 



For the Xe III ion we also included calculations which 
use the basis states of the Xe IV ion (column E{N — 1) 
of Table FVTTT)) . The purpose of these calculations will be 
explained in the negative ions section below. 
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TABLE IV; Ground state removal energy (a.u.), excitation 
energies (cm~^) and g-factors of lowest states of Xe Vll ; 
comparison with experiment 



State 


J 


Expt [26] 


Calculations 










E 


E 


g 


A 









-7.26 




-7.27 




5s5p 







96141 


94889 





1252 






1 


100451 


99394 


1.4846 


1057 






2 


113676 


112598 


1.5000 


1078 


5s5p 


1 po 


1 


143259 


146337 


1.0153 


-3078 


5p^ 


3p 





223673 


224343 





-670 






1 


234685 


235008 


1.5000 


-323 






2 


251853 


252607 


1.3027 


-754 






2 


236100 


237129 


1.1962 


-1029 


5p 


s 





273208 


281328 





-8120 


5s5d 




1 


287772 


291855 


0.5000 


-4083 






2 


288712 


292896 


1.1663 


-4184 






3 


290340 


294591 


1.3333 


-4251 


5s5d 




2 


307542 


317647 


1.0015 


-10105 


5s6s 




1 


354833 


358686 


2.0000 


-3853 


5s6s 







361671 


364853 





-3182 


5p5d 




2 


393792 


398186 


0.7405 


-4394 






3 


401413 


406187 


1.0990 


-4774 






4 


412567 


417863 


1.2500 


-5296 



TABLE V: Ground state removal energy (a.u.), excitation 
energies (cm~^) and (;-factors of lowest states of Xe VI ; com- 
parison with experiment 



State 


J 


Expt [26] 


Calculations 










E 


E 


9 


A 


5s'^5p 




1/2 


-9.71 


-9.72 


0.6667 








3/2 


15599 


15590 


1.3333 


9 


5s5p^ 


4p 


1/2 


92586 


90191 


2.6301 


2395 






3/2 


100378 


97787 


1.7247 


2591 






5/2 


107205 


105036 


1.5629 


2169 


5s5p^ 




3/2 


124870 


125900 


0.8243 


-1030 






5/2 


129230 


129897 


1.2366 


-667 


5s5p^ 


2p 


1/2 


141837 


145429 


1.1878 


-3592 






3/2 


159112 


162903 


1.3119 


-3791 


5s5p^ 




1/2 


157996 


161647 


1.5155 


-3651 


5s^5d 




3/2 


180250 


186188 


0.8058 


-5938 






5/2 


182308 


188093 


1.2004 


-5785 


5s^6s 




1/2 


223478 


224641 


1.9998 


-1163 


5p^ 




3/2 


232586 


232997 


1.3377 


-411 



of the 5s^5p^8s of Xe I can be calculated with good ac- 
curacy with only four states in the basis: 5s, 5pi/2, 5p3/2 
and 8s (see Table |X|) . There are many lower states of 
same parity and total momentum J but we can easily 
get rid of them by not including corresponding single- 
electron states into the basis. 



B. Negative ions 

An interesting question is whether method presented in 
this paper can be used to calculated states of a negative 



TABLE VI: Ground state removal energy (a.u.), excitation 
energies (cm~^) and g-factors of lowest states of Xe V; com- 
parison with experiment 



State 




J 


Expt 126] 


Calculations 










E 


E 


g 


A 


5s'' 5^"^ 


3p 





-11.7 


-11.72 










1 


9292 


8969 


1.5000 


323 






2 


14127 


14643 


1.3744 


-516 


5s^5p^ 




2 


28412 


30169 


1.1256 


-1757 


5s^5p^ 







44470 


47061 





-2591 


5s5p^ 


5go 


2 


92183 


88033 


1.9744 


4150 


5s5p^ 




1 


115286 


115554 


0.6192 


-268 






2 


116097 


116202 


1.2256 


-105 






3 


119919 


120152 


1.3329 


-233 


5s5p^ 


3po 





133408 


134320 





-912 






1 


134575 


135493 


1.4078 


-918 






2 


134703 


135579 


1.3152 


-876 


5s5p^ 




2 


145807 


147030 


1.1261 


-1223 


5s5p^ 


3go 


1 


155518 


160672 


1.7362 


-5154 


5s^5p5d 


3po 


2 


156507 


159419 


0.7036 


-2912 






3 


160630 


163534 


1.0901 


-2904 






4 


169799 


172418 


1.2500 


-2619 


5s5p^ 


Ipo 


1 


169673 


175704 


1.1706 


-6031 



TABLE VII: Ground state removal energy (a.u.), excitation 
energies (cm~^) and g-factors of lowest states of Xe IV 



State 




J 


Expt [26] 


Calculations 










E 


E 


g 


A 


5s''5p=* 




3/2 


-13.2 


-13.27 


1.8987 




5s^5p^ 




3/2 


13267 


14619 


0.9778 


-1352 






5/2 


17511 


18938 


1.2000 


-1427 


5s2 5p^ 


2po 


1/2 


28036 


30149 


0.6667 


-2113 






3/2 


35650 


37446 


1.2569 


-1796 


5s5p'' 


4p 


5/2 


99664 


99466 


1.5814 


198 






3/2 


106923 


106710 


1.7055 


213 






1/2 


109254 


109169 


2.6286 


85 


5s5p'' 




3/2 


121929 


124529 


0.8925 


-2600 






5/2 


125475 


128117 


1.2153 


-2642 


5s^5p^5d 


2p 


3/2 


133027 


135880 


0.8912 


-2853 






1/2 


136796 


139997 


0.7788 


-3201 


5s^5p''5d 


4p 


3/2 


134981 


137617 


0.8501 


-2636 






5/2 


136496 


139103 


1.1064 


-2607 






7/2 


141625 


144013 


1.2649 


-2388 






9/2 


145991 


148958 


1.3088 


-2967 


5s^5p''5d 


2F 


5/2 


141824 


145598 


0.9889 


-3774 






7/2 


145011 


148526 


1.2658 


-3515 


5s^5p^5d 




1/2 


145106 


147933 


0.4037 


-2827 






3/2 


146207 


148762 


1.1487 


-2555 






5/2 


148685 


151840 


1.1799 


-3155 






7/2 


155864 


159785 


1.2361 


-3921 


5s5p'* 


^S 


1/2 


150737 


154437 


1.5219 


-3700 


5s^5p^6s 


4p 


1/2 


157205 


161777 


2.3294 


-4572 






3/2 


165280 


167775 


1.6017 


-2495 






5/2 


170490 


174875 


1.4873 


-4385 


5s25p25d 


4p 


5/2 


159643 


165187 


1.5539 


-5544 






3/2 


161435 


169550 


1.6737 


-8115 






1/2 


162867 


169085 


2.4382 


-6218 
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TABLE VIII: Ground state removal energy (a.u.), excitation 
energies (cm^^) and g- factors of lowest states of Xe III 

State J Expt [26] Calculations 









E 


g 


E 


g 


A 


£;(7V- 1) 






2 


-14.38 




-14.36 


1.4523 




-14.38 









8130 




8313 





-183 


8319 






1 


9794 




9638 


1.5000 


156 


9528 






2 


17099 




19086 


1.0477 


-1987 


17879 









36103 




37280 





-1177 


37392 


5s5p^ 


3po 


2 


98262 




98847 


1.4986 


-585 


98729 






1 


103568 




104334 


1.4553 


-766 


104369 









108334 




108562 





-228 


109016 


5s^5p^5d 


5-^0 


3 


111605 




110836 


1.4702 


769 


113883 






2 


111856 




111066 


1.4623 


790 


114075 






4 


112272 




111366 


1.4813 


906 


114347 






1 


112450 




111506 


1.4882 


944 


114360 









112694 




112142 





552 


114544 


5s^5p^5d 


3-^0 


2 


117240 




118575 


1.8400 


-1335 


122839 






3 


121230 




122482 


1.3114 


-1252 


126365 






1 


121923 




123251 


0.6428 


-1328 


127257 


5s^5p^5d 


Ipo 


1 


119026 




120856 


1.0670 


-1830 


122930 


5s^5p^6s 


5go 


2 


121476 


1.95 


119002 


1.1889 


2474 


125150 


5s^5p^5d 


3po 


2 


124691 




127328 


0.8030 


-2637 


130640 






3 


126120 




129045 


1.0994 


-2925 


132337 






4 


130174 




133970 


1.1386 


-3796 


138195 


5s^5p^6s 


3go 


1 


125617 


1.77 


124152 


1.7102 


1465 


132504 


5s^5p^5d 




4 


127782 




131307 


1.1977 


-3525 


134519 






3 


128349 




132922 


0.8245 


-4573 


137183 






5 


132160 




136422 


1.2000 


-4262 


140435 



TABLE IX: Ground state removal energy (a.u.), excitation 
energies (cm~^) and g-factors of lowest states of Xe II 



State 




J 


Expt [26] 


Calculations 










E 


g 


E 


g 


A 


5s^5p'' 




3/2 


-15.16 




-15.09 


1.3333 








1/2 


10537 




10763 


0.6667 


-226 


5s5p'' 




1/2 


90874 


2.02 


91700 


2.0423 


-826 


5s^5p*6s 


[2] 


5/2 


93068 


1.56 


91729 


1.5639 


1339 






3/2 


95064 


1.38 


94188 


1.4084 


876 


5s^5p''5d 


[2] 


5/2 


95397 


1.36 


95802 


1.3473 


-405 






3/2 


96033 


1.18 


96534 


1.1847 


-501 


5s^5p''5d 


[3] 


7/2 


95438 


1.42 


95783 


1.3940 


-345 






5/2 


106475 




108559 


1.0537 


-2084 


5s^5p''5d 


[1] 


1/2 


96858 


0.50 


97388 


0.5457 


-530 






3/2 


105313 


1.15 


107286 


1.0869 


-1973 


5s^5p''5d 


[4] 


9/2 


99405 


1.31 


100848 


1.3093 


-1443 






7/2 


101536 


1.11 


103581 


1.1524 


-2045 


5s^5p''6s 


[0] 


1/2 


101157 


2.43 


100700 


2.3677 


457 


5s^5p*6s 


[1] 


3/2 


102799 


1.59 


101988 


1.5811 


811 






1/2 


106906 


1.79 


108148 


1.9490 


-1242 


5s^5p''5d 


[1] 


1/2 


104250 


0.56 


104264 


0.6408 


-14 






3/2 


107904 


1.20 


109217 


1.3762 


-1313 


5s^5p''5d 


[0] 


1/2 


105948 


1.36 


107566 


1.1311 


-1618 


5s^5p*6p 


[2]° 


3/2 


111792 


1.61 


112248 


1.6084 


-456 






5/2 


111959 


1.47 


112286 


1.4934 


-327 


5s'^5p*6p 


[3]° 


5/2 


113512 


1.28 


114041 


1.2350 


-529 






7/2 


113705 


1.40 


114266 


1.3984 


-561 


5s^5p''6p 


[1]° 


1/2 


113673 


1.50 


114350 


1.5358 


-677 






3/2 


116783 


1.37 


117621 


1.3609 


-838 



TABLE X: Energies (cm ^) and g-factors of the 5s^5p^8s 
configuration of Xe I 



State 


J 


Expt. 

E 


[26j 

g 


Calculations 
E g 


8s [3/2]° 


2 


90805 


1.465 


92288 


1.5000 




1 


90933 


1.182 


92414 


1.1700 


8s' [1/2]° 









104063 







1 


101426 




104118 


1.3300 



ion. At first glance the answer is no because we use V 
states for the basis and negative ions are not bound in 
the V'^ approximation. However, we may consider the 
following question: what is going to happen if we add one 
more electron to the CI calculations for a neutral atom, 
when basis corresponds to the neutral atom? 

For atoms like xenon, which don't form negative ions, 
it makes more sense to consider more general question: 
can a basis calculated for a system of — 1 valence elec- 
trons be used to calculate many electron states of a sys- 
tem of N electrons? This can be easily checked. Take, 
for example, the Xe IV ion and add one more electron 
in the CI calculations to get the states of Xe III. We've 
done this without adding any new states into the basis. 
Results are presented in the last column of Table IVIIII 
We can see that the results for the ion with the basis of 
the other ion are almost as good as with its own basis. 
Accuracy is a bit lower which is a natural consequence 
of the worsening of the basis. Adding more states to the 
basis would most certainly improve the results. 

This findings are not very surprising since we know 
that any basis set can be used in the CI calculations. 
For example, in Ref. [HI calculations for neutral Kr were 
performed with the basis corresponding to Kr IX! The 
only question is how many states we need to include to 
get reasonable results. It turns out that at least in the 
case of just one more electron there is no need to greatly 
increase the basis. This means that we can also calculate 
states of negative ions by using basis states of a neutral 
atom. 



VIII. CONCLUSION 

In this paper we present a method of calculation for 
many-electron atoms with open shells which is both ac- 
curate and very efficient. The method is based on the 
so called I/^~^^ approximation in which calculations 
start from the highly charged ion with all valence elec- 
trons removed. High accuracy is achieved by inclusion of 
core-valence correlations by means of MBPT. Dominat- 
ing chains of higher order diagrams are included in all 
orders. High efficiency of the method is mostly due to 
the compact basis set for the states of valence elec- 
trons. The method is expected to work well for atoms in 
which valence electrons form a separate shell (defined by 
the principal quantum number) . This is usually the case 
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if valence electrons in the atomic ground state occupy s 
and/or p states. This covers roughly half of the periodic 
table of elements. Calculations for xenon and its ions 
illustrate the use of the method. 
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APPENDIX A: EFFICIENT WAY OF 
CALCULATING E 

The correlation correction operator Si is defined in 
such are way that its overage value over a wave function 
of a valence electron is the correlation correction to the 
energy of this electron: 



where r< = min(r, r') and r> 
need a p function: 



max(r, r'). We will also 



(5e„ = {v\Ysi\v) 



(Al) 



We use the following form of the single-electron wave 
function 



Then the expression (jAip becomes 



Je„ = / / /„(r)E//(r,r')/„(r')drdr' 
fy{r)T.fg (r, r')g^{r')drdr' 
9v{r)^gf{r, r')f^{r')drdr' 
ir)'Egg (r, r')g^ {r')drdr' . 



(A2) 



(A3) 



Note factors and in all terms except the first one. 
These factors make corresponding contributions to be 
very small. Therefore we don't usually include them. 
Only Tiff will be considered in this appendix and we will 
omit the indexes. 



1. Second-order E 

Good efficiency in calculating of S is achieved by di- 
viding the calculations into two steps: 

1. First, all relevant Coulomb Y functions are calcu- 
lated and stored on disk. 

2. Then, S is calculated using stored F-functions. 
The Coulomb F-function is defined as 



' > 

a^9n{r')gm.{r'))dr\ 



Pji{r) = fj{'r)Ii{r) + o? 9j{r)gi{r)). 



(A5) 



Our typical coordinate grid consists of about 1000 points. 
Usually all of them are used to calculate y-functions 
(|A4p . However, there is no need to keep all points for 
the Y and p-functions for consequent calculations. It 
turns out that very little lose of accuracy is caused by 
the use of a subset of points defined as every 4th point 
in the interval 

l/Z <r< Rcore- 

By cutting off the point on short and large distances and 
using only every 4th point in between we reduce the num- 
ber of points by an order of magnitude. Then, Coulomb 
integrals are calculated in an extremely efficient way 

qk{jlmn) = Pji{r{)Ykmn{ri)w^. (A6) 

i=l 

Here p ~ 100 is number of points on the sub-grid and 
Wi are weight coefficients corresponding to a particular 
method of numerical integration. Note that only one 
of two integrations for Coulomb integrals is done on a 
reduced sub-grid. First integration (|A4[) is done with the 
use of all points. 

An expressions for Yj\ via K-functions is 



Si(r,r') = 

/„(r)Yfco„(r)yfc„„(r')/„(r') 

- T,amnk,k2 C2{kik2vamn) X 
/„(r)nio^(r)n,„„(r')/„(r') 

-T,amnk Ciikvabui) X 

Ca+e6-e„— e™ 
+ Eam„fei/c2C4(fclfc2l'a6m) X 

«a+eb-«ii-em 



(A7) 
(A8) 
(A9) 
(AlO) 



(A4) 



Here Ci, C2, C3, C4 are angular coefficients. Expres- 
sions for them can be found elsewhere [IB]- Formulas 
(|A7l) .([A8 |) . (jA9|) . (|A10p . correspond to dia grams 1,2,3,4 on 
Fig. [21 Si is a matrix of size p w 100 in coordinate space. 
Matrix elements of Si are calculated by 

{v\±i\w) = 

Ef=ij=i fv{ri)^iin,rj)fy,{rj)w,Wj. (All) 

Note that we use a two-step procedure to calculate matrix 
elements of Si. First, Si matrix which is independent 
on valence functions is calculated, then matrix elements 
of Si are calculated. To use the same approach for S2 is 
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impractical. As can be seen from Fig. [4] to make E2 in- G is Green function 
dependent on valence states one would have to make ma- 
trices of dimensions 2, 3 and 4. Therefore we just calcu- 
late matrix elements of S2 via Coulomb integrals. Corre- (^1 ~ ^)G{r, r ) = —S{r — r ), 
sponding expressions can be found in Ref . [16|] . Coulomb 

integrals are calculated as in (IA6I) . , ^ . ^ , , . , 

and Q is Coulomb mteraction 



2. All-order S 

We use Feynman diagram technique to include higher- 
order correlations into direct part of Ei (diagrams 1 and 
3 on Fig. [5]) . Corresponding expression is 

/duj 
— G.,j{€ + Uj)Qim^imn{uj)Qnj ■ 

I LI It 

(A12) 

Here 11 is "screened polarization operator" 

n = n[i-gn]-\ 

n is polarization operator 

ea+^) + G{€a - Cj)]l/'a, 

a 



The details of calculation of E^""-' can be found else- 
where 0, Q. Here we only want to mention that op- 
erators H, H, G and Q are matrices of size /i sa 100 in 
coordinate space. Therefore calculation of E^""-* involves 
manipulation of matrices of relatively small size. If we 
also recall that Ei does not depend on valence states and 
needs to be calculated only once then the efficiency of its 
calculation is quiet satisfactory. 

Higher-order correlations in exchange diagrams for Si 
(diagrams 2 and 4 on Fig. ^ and for all diagrams for 
E2 are included via screening factors as explained in the 
text. 
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